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(54) Global positioning system self calibration attitude determination 



(57) An improved GPS Attitude Receiver for space- 
craft uses concurrent line-bias estimates derived from 

GPS signals for attitude determination. Attitude deter- 
mination is obtained from GPS measurements by an in- 
teractive solution of the known "Attitude Master Equa- 
tion", using RF differential phase measurements and 



calibrated parameters of baseline vectors and line-bias- 
es. A "weighted fit error" W, is derived from the sum 
squared of discrepancies between predicted and actual 
measurements, weighted by measurement variances. 
By minimizing W, the line-bias estimations obtained give 
more current data including thermal effects and elimi- 
nate separate ground calibration tests. 
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Description 

[0001] The present application is related to co-pending United States of Annerlca patent application Serial No. 
60/059.092 entitled "GLOBAL POSITIONING SYSTEM SELF CALIBRATION ATTITUDE DETERMINATION", and 
5 based on which priority is herewith claimed and the disclosure of which Is incorporated herein by reference in its entirety. 
[0002] The present invention relates to spacecraft attitude determination systems and more particularly to a system 
that provides pitch, roll, and yaw attitude control for an orbiting satellite using Global Positioning System (GPS) signals 
as Input during operation. 

[0003] A major issue encountered during the development of a suitable GPS Attitude Receiver for spacecraft is the 

10 need for accurate estimates of the electrical delay properties of the system, such as phase lag in the connecting cables 
between the antenna and receiver, which delays are referred to as "line-biases". Existing works relating to the GPS, 
include: 

1) KNIGHT, D., ".A New Method of Instantaneous Ambiguity Resolution", Presented at the Institute of Navigation 
75 7th International Technical Meeting, Salt Lake City, UT Sept. 20-23, 1994; 

2) COHEN, C. "Attitude Determination Using GPS", Ph.D. Dissertation, Stanford University, December 1992; 

3) PARKINSON, B., ET AL, ed., "Global Positioning System: Theory and Applications", Vols. I and II, American 
Institute of Aeronautics and Astronautics, 1996; and 

4) BROCK, J.K., ET AL, "GPS Attitude and Orbit Determination for Space", ION GPS-94 Conference Proceedings, 
20 Salt Lake City, UT, Sept. 20-23, 1 994. 

[0004] These articles have presented some approaches for estimating the line-biases for a stationary platform, that 
offer satisfactory solutions if a suitable test is possible and the line-biases do not change overtime. However, as several 
situations can occur where the test is not possible or the line-biases can be expected to change significantly over time, 
25 such as due to changing thermal conditions, these approaches are somewhat limited. 
[0005] Other articles of interest regarding the GPS in the prior art are: 

5) BROCK, J.K., ET AL, "GPS Attitude Determination and Navigation Flight Experiment", ION GPS-95 Conference 
Proceedings, Palm Springs, CA, Sept. 12-15, 1995; 

30 6) FULLER, R.A.. ET AL. "Spacecraft Guidance and Control with GPS Tensor", Presented at 19th Annual AAS 

Guidance and Control Conference, Breckenridge. CO. Feb. 7-11, 1996; and 7) FULLER, R.A.. ET AL. "GPS Attitude 
Determination From Double Difference Differential Phase Measurements", ION GPS-96 Conference Proceedings, 
Kansas City, MO, Sept. 17-20, 1996. 

35 [0006] It will be seen, however, upon studying all of the foregoing literature, that none involves an approach that 
would produce concurrent generation of attitude parameters, and particularly none using differential phase measure- 
ments obtained from an antenna array while performing phase calibration and reliably estimating attitude system line- 
biases. 

[0007] Accordingly a problem in the art to accurately estimate attitude system line-biases and achieve an approach 
40 that is effective for concurrent generation of attitude using differential phase measurements obtained from an antenna 
array while performing phase calibration. 

[0008] The present invention seeks to provide a method and means for producing concurrent generation of attitude 

using differential phase measurements obtained from an antenna array while performing phase calibration. 

[0009] The invention also seeks to provide a novel approach for accurately estimating attitude system line-biases in 

45 performing the phase, calibration. 

[0010] According to one aspect of the invention there is provided a system for attitude determination in a spacecraft 
using GPS signals to obtain differential phase measurements and using concurrent line-bias estimates, comprising 
antenna means on the spacecraft, including a Master Antenna, for receiving GPS signals, means responsive to said 
received GPS signals, for performing a line of sight, LOS, calculation on the system position parameters, x-y-z, derived 

50 from the received GPS signals, and producing LOS signals indicative thereof, update means, responsive to the LOS 
signals, for performing an attitude update processing and producing an attitude transformation signal, A, baseline 
means for determining the baseline geometry of the locations of the antenna means and producing signals Indicative 
of the baselines, means for producing signals indicative of integer wavelengths in the differential range kij, derived 
from factors comprising initial integer resolutions, line bias estimations, "wraparound" from phase monitoring^ and an 

55 integrity check, means, responsive to GPS signals received by the antenna means, for producing signals indicative of 
differential phase measurements, A^i'i means, responsive to the differential phase measurements, A(J>ij, for producing 
• attitude correction and altitude i-ate correction signals. 68, 50*, 69, SG*, Sy, Sy*. including a line bias correction 6Pj, s'aid 
means comprising means>fdp performing a weighted least squares fit, using ^weighted fit error" W derived from the 
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"Attitude Master Equation', on signals input thereto including, the differential phase measurements signals, A(j)ij, the 
attitude transformation signal. A; the baseline signals, Xj and the signals indicative of integer wavelengths in the dif- 
ferential range, kij, and for producing the line bias correction Apj to provide a concurrent line bias estimate, and means 
for inputting the attitude correction and attitude rate correction signals. 56. 59*. 5o. 5(|)*, 6\|/, 5\i/*. including a line bias 
s correction APj, to the update means for combining with the LOS signals, to obtain an output indicative of attitude and 
attitude rate in pitch, roll, and yaw, 6, 9*, cj), (j)*, \|/*, using concurrent line-bias estimates. 

[0011] According to another aspect of the invention there is provided a method for attitude determination in a space- 
craft using GPS signals to obtain differential phase measurements from which concurrent line-bias estimates are de- 
rived, comprising the steps of: receiving GPS signals with a plurality of antennas on the spacecraft, including a Master 
70 Antenna, producing differential phase measurement signals, A(|)ij. in response to the received GPS signals; 

deriving system position parameters, x-y-z, from the received GPS signals and performing a line of sight, LOS, 
calculation on the position parameters, x-y-z, and producing LOS signals indicative thereof; in response to the 
LOS signals, performing an attitude update and producing an attitude transformation signal, A, determining the 

75 baseline geometry of the locations of the antennas and producing signals indicative of the baselines, Xj; 

producing signals indicative of integer wavelengths in the differential range, kij, derived from factors comprising 
initial integer resolutions, line bias estimations, "wraparound" from phase monitoring, and an integrity check; 
performing a weighted least squares fit, using "weighted fit error" W derived from the "Attitude Master Equation", 
on the differential phase measurement signals, A<l)ij, the attitude transformation signal. A, the baseline signals, Xj. 

20 and said signals indicative of integer wavelengths in the differential range, kij, and producing a line bias correction 

APj to provide a concurrent line bias estimate; 

applying the line bias correction A^j to the differential phase measurement signals, A(|)ij, and producing attitude 
correction and attitude rate correction signals, 69, 59*. 6(t> 5ct>*, 6\|/, Sy*, including the line bias correction A^j; 
combining the LOS signals with the attitude correction and attitude rate correction signals, 59, 59*, 59, 5<|)*, 5\j/, 
25 S\j/*^ and producing an output indicative of attitude and attitude rate in pitch, roll, and yaw, 9, 9*. (j), 0*, y*, with 

concurrent line-bias estimates. 

[0012] According to another aspect of the invention there is provided a GPS attitude receiver in a spacecraft for 
determining spacecraft attitude using GPS signals to obtain differential phase measurements and using concurrent 

30 line-bias estimates, comprising: means on the spacecraft for receiving GPS signals from a plurality of spacecraft an- 
tennas, including a Master Antenna, for receiving GPS signals, means for deriving spacecraft position parameters, x- 
y-z, from the received GPS signals; means for performing a line of sight, LOS, calculation on the spacecraft position 
parameters, x-y-z, and producing LOS signals indicative thereof, update means, responsive to the LOS signals, tor 
performing an attitude update processing and producing an attitude transformation signal. A, baseline means for de- 

35 termining the baseline geometry of the locations of said plurality of antennas and producing signals indicative of the 
baselines, Xj, means for producing signals indicative of integer wavelengths in the differential range^ kij, derived from 
factors comprising initial integer resolutions, line bias estimations, "wraparound" from phase monitoring, and an integrity 
check, means, responsive to the GPS signals received, for producing signals indicative of differential phase measure- 
ments, A(t)ij, means, responsive to the differential phase measurements, Ac^ij, for producing attitude correction and 

40 attitude rate correction line signals, 59, 59*, 59, 5(t)*, Sy, 5vi/*, including a line bias correction A^j, the means comprising, 
means for performing a weighted least squares fit, using "weighted fit error" W derived from the "Attitude Master Equa- 
tion", on signals input theretO: including the differential phase measurements signals. Aoij, the attitude transformation 
signal A, the baseline signals,, Xj, and the signals indicative of integer wavelengths In the differential range kij, and for 
producing the line bias correction Apj to provide a concurrent line bias estimate, and means for inputting the attitude 

45 correction andattitide rate correction signals. 59, 59*. 5(J), 5(t)*, 6\j/, 5\i/*, including a line bias correction Apj, to the update 
means for combining with the LOS signals, to obtain an output indicative of attitude and attitude rate in pitch, roll, and 
yaw, 9, e*. 0*. Y, xi/*, using concurrent line bias estimates. 

[0013] The invention enables the provision of a method and means for significantly extending the usefulness of 
systems that accomplish attitude determination using signals from the Global Positioning System (GPS), such as a 

50 system using the technique described in the previously-noted paper 1 ) by D. KNIGHT, "A New Method of Instantaneous 
Ambiguity Resolution", and particularly the invention enables the provision of an improved GPS Attitude Receiver for 
spacecraft. In contrast to the prior art pre-calib ration, fixed parameter, "'Self -Survey" technique developed by G. COHEN 
and described In the above-noted technical paper 2), attitude determination from GPS measurements are obtained by 
an interactive solution of the known "Attitude Master Equation", set forth, using RF differential phase measurements 

55 and calibrated parameters of baseline vectors and line-biases. Baselines can be determined mechanically, but line- 
biases require special calibration, and line-biases in orbital operation are utilized on a nearly continuous basis. Con- 
sequently, estimating line-biases accurately in a spacecraft GPS Attitude receiver poses a considerable challenge, 
aiid the present Invention permits generation of concurrent accurate estimates suitable for this purpose. The estlma- 
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tions obtained with the invention give more current data including thernnal effects and eliminate separate ground cali- 
bration tests. 

In obtaining the estimates, the phase measurements and calibrated parameters of baseline vectors and line-biases to 
be used are related In the "Attitude Master Equation" as follows: 



Ac|)ij = Si'^A'^Xj-klj+pj + Vij. 



where: 

10 

A(t)ij = the Differential Phase Measurement with respect to the Reference or Master Antenna signal; 

i = GPS Channel, 1, 2, 9; 

j = Antenna Baselines, 1 , 2, 3; 

S'\^ = the GPS Satellite Line of Sight (LOS) vector In Earth Fixed Earth Centered, ECEF, coordinates; 

15 AT = the Transformation Matrix, ECEF from Body Coordinates; 

Xj = the Antenna Baseline Vector In Body Coordinates; 

klj = the Differential Integer Ambiguity; 

Pj = the Differential Line-Bias or electrical path length; and 

Vi] = the Phase Measurement Noise. 

20 

[0014] The Differential Range Measurement is the distance in the direction of the wavefront of the slave from the 
Reference or Master Antenna, i.e.,: 



25 Arij = Si^ A^ Xj = A^ij + klj - Pj. 

[001 5] An alternate state space formulation of this Attitude Master Equation, as set forth In the above-noted technical 
paper by D. KNIGHT is: 

30 ei = Hix-nl + el, 1= 1...m 

e = Hx-n +^ 

X = [a / b] or transposela bl, H = ILOS1 1 1 

35 

where: 

e = Carrier phase measurement; 

H = Observation matrix; 
40 X = State, typically antenna position and clock bias; 

n = the integer combination, one per satellite and there are m satellites, numbered i = 1 ...m, with each ni = from 

1...to m representing a carrier cycle Integer and associated uncertainty bounds; 

Gi = phase measurement error with an assumed Gaussian probability distribution and variance Rl; 

a = antenna position; 
45 b = clock bias or electrical path lengthy; and 

LOS = llne-of-sight, the direction of incident signal arrival. 

[0016] All units of distance, phase, and time are In carrier wavelengths. 

[0017] Viewing X and n as parameters to be adjusted, the measured phases, 0i, are then independent random var- 
so iables, and the joint probability density of 9 Is the product of the Individual densities. Evaluating the joint probability 

density of 0: 

p(e ) = p(0i ) p(e2) p(e3) ... p(Gm) 

55 * p(ei) = (1/(27cRi)^.5) exp -[(91 - HI x + nl)'^2 / 2Ri] 

[0018] The log-likelihood function Is: ' ife, . 
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In p(e) = X (i=1 to m) - 1/2 ln{27iRi) - 1/2 (Gi - Hi x + ni)^2 / Ri 
and p(9) is maximized when, 

5 [001 9] W = Z (i=:1 to m) (0i - Hi x + ni)^2 / Ri is minimized. W is the "weighted fit error", the sum squared of discrepancies 
between predicted and actual measurements, weighted by measurement variances. The maximum-likelihood estimate 
is the (X, ri) that minimizes W, and thus is the one sought as the most accurate in evaluatng the line-biases. 
[0020] Using the relationships expressed in the foregoing equations with the GPS and other measurement data 
obtainable with the invention, improved attitude determination can be achieved. 

10 [0021] In order that the invention and its various other preferred features may be understood more easily, some 
embodiments thereof will now be described, by way of example only, with reference to the drawings, in which:- 

Figure 1 is a flow diagram illustrating a navigation solution flow in a GPS Attitude Receiver for spacecraft resulting 
in X, y, z, and b parameters for processing in accordance with the present invention, 
^5 Figure 2 is an illustrative diagram of a GPS attitude processing system algorithm in accordance with the invention, 

Figure 3 illustrates an attitude algorithm processing archrtecture for implementing the system of Figure 2. 

The present invention is directed to significantly extending the usefulness of systems for accomplishing attitude deter- 
mination based on signals from the Global Positioning System (GPS) and embodies a method and means for this 
purpose which use and improve upon an approach similar to that described in the above-noted technical paper 1) of 
D. KNIGHT. - 

[0022] "A New Method of Instantaneous Ambiguity Resolution", and provide a GPS Attitude Receiver achieving 
accurate estimates of line-biases, as will be understood from the following detailed description. 
[0023] To begin with, an example of a current advanced aerospace qualified GPS receiver in which preferably the 
present invention may be implemented is the "GPS TENSOR"r"" (TM of the applicant of the present application), such 
as disclosed in the above-noted article 6) by present co-inventor R. A. FULLER, and others, entitled "Spacecraft Guid- 
ance and Control with GPS Tensor", which article is incorporated herein by reference. Updated information thereon is 
disclosed in another article by the present co-inventors, and others, entitled "GPS Tensor, An Attitude and Orbit De- 
termination System for space", R. FULLER, D. HONG. S. HUR-DIAZ, J. RODDEN, M. TSE. presented at the Institute 

^0 of Navigation 10th International Technical Meeting, ION GPS-97, September 16-19, 1997, Kansas City, MO, which 
article is also herein incorporated by reference. The GPS TENSOR^*^ generates navigation position/velocity data, 
precision timing signals, and 3 axes of spacecraft attitude, utilizing 9 channels, and is particularly suited for use with 
the "GLOBALSTAR^"^" (TM of Globalstar, Inc.) communication satellite constellation presently in the process of being 
implemented. The operation of this receiver starts with a coded RF signal from the sum of the present 25 GPS Satellite 

55 Vehicles (SVs) in a constellation in orbit at about 20.000 km, half synchronous altitude. The first step in processing is 
to measure "pseudoranges" between the spacecraft or satellite carrying the receiver and the GPS SVs. The pseudor- 
ange is the time delay or "distance" between the time a signal is sent from a GPS SV and the time of its receipt by the 
spacecraft receiver plus the clock bias offset of the receiver clock, i.e., the actual transmission time, which is indicative 
of true geometrical distance between the SV and receiver at the time of signal transmission plus the effect of the 
receiver clock bias. 

[0024] The pitch, roll, and yaw attitude, 0, (j>, v}/, determination approach of the invention, follows from and improves 
upon that first developed by the previously-mentioned Dr. C. COHEN of Stanford University", and disclosed in technical 
paper 2). In contrast to the fixed calibration of COHEN, with the newly described arrangement navigation calculations 
are available and in progress. These attitude calculations use signals from an SV received by multiple antennas on 
the spacecraft, although in actuality the navigation function only requires one antenna. In fact, a Master Antenna is 
designated from among the multiple antennas for pseudorange determination, and re<iesignation is provided for 
switching to another antenna upon changes in the signal to noise ratios (SNRs) on the tracking channels. The GPS 
TENSORS", of the preferred embodiment, nominally uses signals from 4 rectangularly arrayed antennas to determine 
attitude. The signals are RF signals received at all of the antennas from each SV and from which signals differential 

5^ phase measurements A4)ij are derived. Several SV signals are required. The slight phasing differences of signals ar- 
riving at different times from different sources at the different antennas provide the basis for determining the attitude 
of the spacecraft with respect to the local reference frame. The attitude calculations presently used need some data 
base parameters called "Self Survey Data". These include the vector distances between the antennas and the desig- 
nated Navigation Master Antenna, called "baselines", and the electrical phase losses or "line-biases" in the cable 

^5 lengths from the antennas to the receiver. The original attitude system used to implement the approach first developed 
by Dr. Clark COHEN, noted above, previously mentioned required a special measurement process for determining 
these parameters that took several hours with the fully configured (antenna - cable - receiver) system collecting RF 
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data from GPS SVs (see Dr. COHEN'S U.S.Pat. No.5,548, 293. issued 20 August 1996, for "System and Method for 
Generating Attitude Determinations using GPS") This system was originally applied in aircraft which could be parked 
on a runway for the required measurement periods. However, the application of this Self Survey Data collection method 
in the production of the GLOBALSTAR™ satellites was prohibitively difficult since, although the vector baselines could 
s be found from the mechanical positions of the antennas, no off-line measurements of each of the line-biases were 
available. 

[0025] In addition to the difficulty of predetermining the line-biases, the expectation of significant temperature vari- 
ability of the line-biases in orbiting GLOBALSTAR™ satellites motivated the development of a system and software 
that could estimate the line-biases concurrently with the flight operation. This development has led to the present 
10 invention as will now be described. 

[0026] Attitude determination from GPS measurements generally may be obtained by an interactive solution of the 
"Attitude Master Equation" with RF differential phase measurements and calibrated parameters of baseline vectors 
and line-biases. In obtaining the estimates, the phase measurements and calibrated parameters of baseline vectors 
and line-biases to be used are related in the "Attitude Master Equation" as follows: 

IS 

A(t>ij = Si*"" A*'' Xj - kij + Pj + Vij, 



where: 

20 

A(|)ij = the Differential Phase Measurement with respect to the Reference or Master Antenna signal; 

I = GPS Channel, 1 , 2, 9; 

j = Antenna Baselines, 1 , 2, 3; 

SF = the GPS Satellite Line of Sight (LOS) vector in Earth Fixed Earth Centered, ECEF, coordinates; 

25 AT = the Transformation Matrix, ECEF from Body Coordinates; 

Xj = the Antenna Baseline Vector in Body Coordinates; 

kij = the Differential Integer Ambiguity; 

Pj = the Differential Line-Bias or electrical path length; and 

Vij = the Phase Measurement Noise. 

30 

[0027] The Differential Range Measurement is the distance in the direction of the wavefront of the slave from the 
Master Antenna, i.e.,: 



35 Arij = Si^ A^ Xj = A(|)ij + kij - Pj. 

[0028] While the baselines can be readily determined mechanically line-biases require special calibration since line- 
biases In orbital operation are on a nearly continuous basis. It is this calibration which the present invention is designed 
to generate. The estimation obtained with the invention gives more current data with thermal effects and eliminates 
40 separate ground calibration tests. 

[0029] More particularly, an alternate state space formulation of this Attitude Master Equation, as set forth in the 
above-noted technical paper 1) by D. KNIGHT, is as follows: 



45 



Bi = Hi X - ni + Gi, i = 1...m 
e = H X - n +^ 

X = [a / b] or Iransposela bl, H = ILOSi 11 
so where: 

0 = Carrier phase measurement; 
H = Observation matrix; 

X = State, typically antenna position and clock bias; 
55 n = the Integer combination, one per satellite and there are m satellites, numbered i = 1 ...m. with each ni = from 

1...to m representing a carrier cycle integer and associated uncertainty bounds; 
Gi = phase measurement- error with an assumed Gaussian probability distribution and variance Ri; 
a = antenna position; 
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b = clock bias or electrical path length; and 

LOS = line-of-sight, the direction of incident signal arrival. All units of distance, phase, and tinne are In carrier 
wavelengths. 

5 [0030] Viewing X and n as parameters to be adjusted, the measured phases, ei, are then independent random var- 
iables, and the joint probability density of 0 is the product of the individual densities. Evaluating the joint probability 
density of 9: 

p(e) = p(ei)p(e2) p(e3) ...p(0m) 

10 

p{Gi) = (1/(27rRi)'^.5)exp-[(ei- Hi x + nl)'^2 / 2Ri] 
[0031] The log-likelihood function is; 

In p(9) = £ (i=1 to m) - 1/2 ln(27cRi) - 1/2 (91 - Hi x + ni)^2 / Ri 
and p(e) is maximized when, 

[0032] W = Z (1=1 to m) (0i - Hi X + ni)'^2/ Ri Is minimized. W is the ^weighted fit error", the sum squared of discrepancies 
2^ between predicted and actual measurements, weighted by measurement variances. The maximum-likelihood estimate 
is the (x, n) that minimizes W, and thus is the one sought as the most accurate in evaluating the line-biases. This 
problem formulation is presented in, D. KNIGHT U.S. Patent No. 5,296,861, issued March 22, 1994. 
[0033] The preceding problem statement has been used previously for the determination of the integer number of 
cycles, ni, with known electrical path length measurements. The process implies a dedicated measurement procedure 
for external determination on the electrical paths. The present invention enables expansion of the estimation process 
to include a self determination of the electrical path length which eliminates the need for external measurement pro- 
cedures. In addition, the "real time" estimation automatically compensates for slowly varying electrical path lengths as 
might occur on a thermally varying spacecraft environment. 

[0034] In applying the equation information in the GPS TENSOR^m it should be understood that there are two nav- 

30 igation solutions performed in the GPS TENSOR, i.e., the Standard Position Sen/ice (SPS) and the Navigation Filter 
(NF), as illustrated in the flow diagram of Figure 1 . As seen in the Figure, a Delay Lock Loop (DLL) 10 is provided to 
receive from an SV coded RF signals and measure fractional millisecond parts of the pseudoranges with the integer 
number of milliseconds in a pseudorange being initially calculated upon SV signal acquisition. Both the pseudorange 
and clock bias values are divided into terms comprising Integer millisecond distances, e. g., 300 km. plus a fractional 

55 millisecond segment. Additional integers are rolled over, from the raw Code Phase millisecond measurement, each 
time the fractional millisecond exceeds an integer, and are input at 11. The integer clock bias is similarly reset and 
rolled over at 13 from the bias calculated in either the Standard Position Service (SPS) 12 or the Navigation Filter (NF) 
14 and this clock bias is used at 11 to offset the pseudorange measurement data that goes into both-the SPS and NF 
calculations through the consistency check at 15 which passes the millisecond integer corrected pseudorange. The 

40 MF 14 computations have an independently determined bias estimate so that the HF bias value is incremented at the 
same time the adjustment is made to the SPS 12 computations to keep its clock bias consistent with the adjusted 
pseudorange measurements and vice-versa. The GPS time reference is the standard throughout the network lor com- 
munications systems such as the GLOBALSTAR satellites. The SPS 12 numerically solves for the 4 unknowns, i.e., 
the x-y>z positions and the clock bias b, using 4 pseudorange measurements. This Service (SPS) 12 is subject to 

^5 position calculation outages and without these 4 measurements no SPS solution is possible. The alternate navigation 
solution involving the Navigation Filter 14 algorithm uses individual pseudorange measurements to adjust a continuous 
running dynamic model of the 4 system position parameters, x-y-z and clock bias b, and their derivatives. The Navi- 
gation Filter 14 can use up to all 9 measurement channels available In the GPS TENSOR. It also will generate an 
estimate of the position parameters, even in the event of short term absences of new pseudorange inputs, based on 

5^ dynamic extrapolation of the preceding operation. The NF 1 4 operation requires an initialization period o1 several orbits 
. for its parameters to converge to the proper solution, e.g., 120 SPS solutions, and it estimates 9 states, i.e., 3 position, 
3 velocity, and clock bias, bias rate of change, and bias acceleration. Once the parameters are converged, the NF 
solution is more accurate than the SPS solutkjn since it continually averages more data. The NF solution is also more 
robust since it doesnt suffer outages when there are less than 4 SVs available for processing. The Interaction of these 

55 two solutions in conditioning the pseudorange input data may be dealt with by the melding of the integer millisecond 
roll-over from the Code Delay Lock Loop (DLL) 1 0 with the millisecond resets and the Pulse Per Second (PPS) output 
generated at 1G from the clock bias part 17 of the navigation solution. The PPS output at 16 is more robust with the 
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NF solution than with the SPS solution since, as noted, the latter is more susceptible to outages. For satellite systems, 
such as GLOBALSTAR, the precision timing signal of the PPS output 16 is critical. There are separate output pins for 
these pulses which are produced at a 1 second Interval accurate to a microsecond. The entire network of GLOBALSTAR 
communication frequencies are locked onto the "universal" GPS Time as detected and regenerated by the GPS TEN- 
SOR The generation of the PPS is based on the clock bias parameter of the navigation calculations. The GPS TENSOR 
operates such that the default selection of the NF calculated clock bias is used for PPS generation. The system position 
parameters x-y-z are used in a GPS attitude processing algorfthm. 

[0035] An illustrative diagram of a suitable GPS attitude processing algorithm employed In the invention Is shown in 
Figure 2. Multiple phase measurements, derived from multiple wavelength signals received from a number of available 
GPS SVs, are input at a rate of 10 Hz from the antennas, in the form of differential phase measurements A<j)ij and rate 
measurements A4)*ij, and subjected to a weighted least squares fit at a node 20, using the "weighted fit error" W. Node 
20 also receives an input 21 , indicative of attitude transformation A after an attitude update processing at node 22, a 
baselines Xj input 23 after a baseline determination at node 24, and an integer wavelengths in the differential range 
kij Input 25, which wavelenghts are derived at node 26 from factors such as an Initial Integer resolution, line bias 
estimation, "wraparound" from phase monitoring, and an integrity check. The least squares fit at 20 outputs an indication 
27 of attitude correction and rate, e.g., for pitch, 56, 60', which Includes a line bias correction Apj. Indication 27 is Input 
for attitude update processing at 22 along with an input 29 representing an LOS calculation at node 28 using the x-y- 
z parameters from the solution flow In Figure 1 . which results In the desired pitch, roll, and yaw attitude and rate 
determinations, 6, 9*, <p.. <^*, y, y*. 

[0036] An appropriate attitude algorithm processing architecture for implementing the present invention is shown in 
Figure 3. The received RF GPS signals are passed through a hard limlter 30 and converted or coded into IF signals 
which are input, e.g., to 9 channels, to produce the differential phase and rate measurements A())lj and A^*\] Input to 
node 20 In Figure 2. These A<^\\ and A(|>*ij measurements are input at 31 to an attitude solution unit 32, essentially 
functioning as nodes 20 and 22, which also receives an input 33 indicative of baseline vectors and line biases, and an 
input 35 indicative of integer corrections. The attitude solution unit 32 provides feedback outputs to modify the inputs 
33 and 35, i.e., one output 37 to a line bias estimator unit 34 having an output that is added to the baseline vector and 
line bias input 33; another output 39 to an integer estimator unit 36, and a solutions residuals output 41 to an integrity 
checking unit 38. The output 43 of the integrity checking unit 38 is supplied, along with the output 45 from the integer 
estimator unit 36, to an adder 40 that provides the integer corrections Input 35 to attitude solutions unit 32, In response 
to the inputs at 31 . 33, and 35, the attitude solution unit 32 provides an output 46 indicative of attitude and attitude 
rate, G, 0*. (j), <^*, \|/, for pitch, roll, and yaw in keeping with the algorithm of Figure 2. 

[0037] It will accordingly be seen that a method and means have been described that significantly extends the use- 
fulness of systems for accomplishing attitude determination based on the use of signals from the Global Positioning 
System (GPS). A system so adapted for attitude determination in a spacecraft is dependent only on electrical phase 
measurements coupled with input geometric input data. No separate prior measurement "Self Survey" procedure is 
required for determination of electrical path phase length. Electrical phase lengths need not be constant during oper- 
ation. 



40 Claims 

1. A system for attitude determination in a spacecraft using GPS signals to obtain differential phase measurements 
and using concurrent line-bias estimates, comprising antenna means on the spacecraft, including a Master An- 
tenna, for receiving GPS signals; means (28), responsive to said received GPS signals, for performing a line of 

45 sight, LOS, calculation on the system position parameters, x-y-z, derived from the received GPS signals, and 

producing LOS signals indicative thereof, update means (22), responsive to the LOS signals, for performing an 
attitude update processing and producing an attitude transformation signal, (A), baseline means (24) for determin- 
ing the baseline geometry of the locations of the antenna means and producing signals (Xj) indicative of the base- 
lines, means (26) for producing signals indicative of Integer wavelengths in the differential range kij, derived from 

^0 factors comphsing initial integer resolutions, line bias estimations, "wraparound" from phase monitoring, and an 

integrity check, means, responsive to GPS signals received by the antenna means, for producing signals indicative 
of differential phase measurements, A<|>ij, means (32), responsive to the differential phase measurements, A<|)ij, for 
producing attitude correction and attitude rate correction signals, 50, 60*, 5<|>, 6<|)*, 5v|/, S\|/*, including a line bias 
correction Spj, said means comprising means (20) for performing a weighted least squares fit, using "weighted fit 
error" W derived from the "Attitude Master Equation", on signals input thereto including, the differential phase 
measurements signals, A(|>ij, the attitude transformation signal, A; the baseline signals, Xj and the signals indicative 
of integer wavelengths in the differential range, kij.^S'nd for producing the line bias correction APj to provide a 
concurrent line bias estimate, and means for inputting the attitude correction and attitude rate correction signals, 
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• *> 

89. 59*, 5<|). 5<t)*, Sy, 5v|/*, including a line bias correction APj, to the update means (22) for combining with the LOS 
signals, to obtain an output indicative of attitude and attitude rate in pitch, roll, and yaw, 9, 9*, <t». (^*, \|/*, using 
concurrent line-bias estimates. 

2. A system as claimed in claim 1 . comprising a plurality of tracking channels for receiving the GPS signals and 
means (14) for selecting the Master antenna from among the antenna means based on the signal to noise ratios 
(SNRs) on the tracking channels. 

3. A syistem as in claim 1 or 2, comprising means for minimizing the "weighted fit error" W. 

4. A method for attitude determination in a spacecraft using GPS signals to obtain differential phase measurements 
from which concurrent line-bias estimates are derived, comprising the steps of: receiving GPS signals with a plu- 
rality of antennas on the spacecraft, including a Master Antenna; producing differential phase measurement sig- 
nals, A(|>ij, in response to the received GPS signals; deriving system position parameters, x-y-z, from the received 

is GPS signals and performing a line of sight, LOS, calculation on the position parameters, x-y-z, and producing LOS 

signals indicative thereof; in response to the LOS signals, performing an attitude update and producing an attitude 
transformation signal, A; determining the baseline geometry of the locations of the antennas and producing signals 
indicative of the baselines, Xj; producing signals indicative of integer wavelengths in the differential range, kij, 
derived from factors comprising initial integer resolutions, line bias estimations, "wraparound" from phase moni- 
20 toring, and. an integrity check; performing a weighted least squares fit, using "weighted fit error" W derived from 

the "Attitude Master Equation", on the differential phase measurement signals, A(|)ij, the attitude transformation 
signal, A, the baseline signals, Xj, and said signals indicative of integer wavelengths in the differential range, kij, 
and producing a line bias correction A^j to provide a concurrent line bias estimate; applying the line bias correction 
APj to the differential phase measurement signals, A(i»ij, and producing attitude correction and attitude rate correc- 
ts tion signals, 59, 59*, 5(|>*, 6^, 5vj/*, including the line bias correction APj; combining the LOS signals with the 
attitude correction and attitude rate correction signals, 69, 59*, 5<t>, Scj)*, 5\|/, 5\|/*, and producing an output indicative 
of attitude and attitude rate in pitch, roll, and yaw, 9, 9*, (j), <!>*, \(/, y*, with concurrent line-bias estimates. 

5. A method as claimed in claim 4, wherein the GPS signals are received on a number of tracking channels and the 
30 Master Antenna is designated from among the plurality of antennas based on the signal to noise ratios (SNRs) on 

the tracking channels. 

6. A method as claimed in claim 4 or 5, wherein the "weighted fit error" W is minimized. 

3S 7. A GPS attitude receiver in a spacecraft for determining spacecraft attitude using GPS signals to obtain differential 
phase measurements and using concurrent line-bias estimates, comprising: means on the spacecraft for receiving 
GPS signals from a plurality of spacecraft antennas, Including a Master Antenna, for receiving GPS signals, means 
for deriving spacecraft position parameters, x-y-z, from the received GPS signals; means (28) for performing a 
line of sight, LOS, calculation on the spacecraft position parameters, x-y-z, and producing LOS signals indicative 

40 thereof, update means (22), responsive to the LOS signals, for performing an attitude update processing and 

producing an attitude transformation signal, A, baseline means (24) for determining the baseline geometry of the 
locations of said plurality of antennas and producing signals indicative of the baselines, Xj, means (26) for producing 
signals indicative of integer wavelengths in the differential range, kij, derived from factors comprising initial integer 
resolutions, line bias estimations, "wraparound" from phase monitoring, and an integrity check, means, responsive 

45 to the GPS signals received, for producing signals indicative of differential phase measurements, A(|)!j, means (32), 

responsive to the differential phase measurements, A<l)ij, for producing attitude correction and attitude rate correc- 
tion signals, 59, 59*. 54>, 5(J)*, Svj;, 5\}/*, including a line bias correction APj, the means comprising, means (20) for 
performing a weighted least squares fit, using "weighted fit error" W derived from the "Attitude Master Equation", 
on signals input thereto, including the differential phase measurements signals, A(|jij, the attitude transformation 

so signal A, the baseline signals, Xj, and the signals indicative of integer wavelengths in the differential range kij, and 

for producing the line bias correction APj to provide a concurrent line bias estimate, and means for inputting the 
attitude correction and attitude rate correction signals, 50, 59*, 5<^), 6(|)*, 5\]/, 5\}/*. including a line bias correction APj, 
to the update means (22) for combining with the LOS signals, to obtain an output indicative of attitude and attitude 
rate in pitch, roll, and yaw, 9, 9*, <|), <()*, y, y*, using concurrent fine bias estimates. 



55 



A receiver as claimed in claim 7, further comprising a plurality of tracking channels for receiving said GPS signals 
and means (14) for selecting the Master Antenna froni among the antenna means based on the signal to noise 
ratios (SNRs) on the tracking channels. 
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9. A receiver as claimed in claim 7 or 8, comprising means for minimizing said "weighted fit error" W, 
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(57) An improved GPS Attitude Receiver for space- 
craft uses concurrent line-bias estimates derived from 
GPS signals for attitude determination. Attitude deter- 
mination is obtained from GPS measurements by an in- 
teractive solution of the known "Attitude Master Equa- 
tion", using RF differential phase measurements and 



calibrated parameters of baseline vectors and line-bias- 
es. A "weighted fit error" W, is derived from the sum 
squared of discrepancies between predicted and actual 
measurements, weighted by measurement variances. 
By minimizing W: the line-bias estimations obtained give 
more current data including thermal effects and elimi- 
nate separate ground calibration tests. 
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